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We investigate the classical stability of higher-dimensional Schwarzschild black holes with 
respect to linear perturbations in the framework of a gauge-invariant formalism for gravi- 
tational perturbations of maximally symmetric black holes. This formalism was recently 
developed by the present authors. The perturbations are classified into three types, those 
of tensor, vector, and scalar modes, according to their tensorial behaviour on the spherical 
section of the background metric. The vector- and scalar-type modes correspond, respec- 
tively, to the axial and polar modes in the four-dimensional case. We show that for each 



in 



mode of the perturbations, the spatial derivative part of the master equation is a positive, 

■ self-adjoint operator in the L 2 -Hilbert space, and hence that the master equation for each 
type of perturbation has no normalisable negative modes that would correspond to unstable 
solutions. 

■ In the same Schwarzschild background, we also analyse static perturbations of the scalar 
£f~^ ' mode and show that there exists no static perturbation that is regular everywhere outside the 

event horizon and is well behaved at the spatial infinity. This confirms the uniqueness of the 
, higher-dimensional spherically symmetric, static vacuum black hole, within the perturbation 

i-Q ■ framework. 

Our strategy for treating the stability problem is also applicable to other higher- 
CIh, dimensional maximally symmetric black holes with a non- vanishing cosmological constant. 

We show that all possible types of maximally symmetric black holes (including the higher- 
• dimensional Schwarzschild-de Sitter and Schwarzschild-anti-de Sitter black holes) are stable 

with respect to tensor and vector perturbations. 

X 

5_i ■ §1. Introduction 

Higher-dimensional theories of gravity have recently attracted much attention, 
and black hole solutions have played an important role in revealing various as- 
pects of gravity peculiar to higher-dimensional systemsP^^® A noteworthy 
fact is that, compared to the severely constrained situation in four spacetime di- 
mensions P'lHJQ higher-dimensional black holes can have many very different va- 
rieties For this reason, it is important to study the dynamics of 

fieidPunnianiEj and metric perturbationf pfl3m 

in such black hole back- 
grounds and to find stable solutions. Amongst the various types of black hole solu- 
tions, a natural higher-dimensional generalisation of the Schwarzschild metric — also 
known as the Schwarzschild- Tangherlini metric^ — has been assumed to be (or be- 
lieved to be) stable, like its four-dimensional counterpart, and for this reason it has 

*' E-mail: A.Ishibashi@damtp.cam.ac.uk 
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been conjectured that this metric describes a possible final equilibrium state of clas- 
sical evolution of an isolated system, or a certain class of unstable black holes in 
higher dimensions. However, the question of the stability of a Schwarzschild black 
hole itself is not a trivial matter, and, to our knowledge, it has not yet been fully 
analysed in higher dimensions. 

In the four-dimensional case, the perturbation analysis of the Schwarzschild 
metric was first carried out by Regge and Wheeler and its stability was essentially 
confirmed by Vishveshwarap^ followed by further detailed studies The basic 

observation of those works is that unstable perturbations, whose growth in time 
is unbounded, are initially divergent at the event horizon (provided that they are 
vanishing at large distances) and therefore are physically unacceptable. 

A key step in the analysis is the reduction of the perturbation equations to a 
simple, tractable form. For this purpose, in the case of a four-dimensional spacetime, 
the spherical symmetry of the background Schwarzschild metric is utilised to classify 
the metric perturbations into two types, axial (odd) modes and polar (even) modes, 
according to their behaviour under the parity transformation on the two-sphere. 
To linear order, the perturbation equations for these two types of perturbations 
are independent, and therefore they reduce to a set of two Schrodinger-type 2nd- 
order ordinary differential equations (ODEs), the Regge- Wheeler equation for axial 
perturbations and the Zerilli equation^ «^5I f or polar perturbations. 

In the higher-dimensional case, the role of the two-sphere is played by an n- 
sphere: According to their tensorial behaviour on the n-sphere, there are three types 
of perturbations, tensor, vector, and scalar perturbations. The first of these is a 
novel type of perturbation that exists only in the higher-dimensional case, while 
the last two types correspond, respectively, to the axial and polar modes in the 
four-dimensional case. 

In order to verify the stability of higher-dimensional Schwarzschild black holes, 
one must investigate their stability with respect to all the types of perturbations 
mentioned above. However, as far as the present authors know, all analysis made 
to this time has been carried out for only tensor perturbations pi and the situation 
for vector and scalar perturbations has not yet been investigated: The stability of 
Schwarzschild black holes in higher dimensions is not yet established. The treatment 
of scalar perturbations*) is the most difficult part of this analysis, because the lin- 
earised Einstein equations for a scalar perturbation are much more complicated than 
those for a tensor perturbation and do not reduce to the form of a single 2nd-order 
Schrodinger-type ODE in the higher-dimensional case. 

Recently, this reduction was carried out by the present 

authorPJ in the back- 
ground of maximally symmetric black holes with a non-vanishing cosmological con- 



*' The term "scalar perturbation" is often used in reference to a perturbation of a free scalar 
field in works involving quasi-normal mode analysis of black holes. This is conceptually quite 
different from the scalar perturbations considered in the present paper, which are modes of metric 
perturbations of a vacuum spacetime. Furthermore, the effective potentials for these two cases are 
different. Hence, although the quasi-normal mode analysis of a test scalar field^^^^'^' may 
shed some light on the stability problem of black holes in higher dimensions, an exact analysis of 
the problem based on the formulation given in the present paper is necessary. 
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stant in an arbitrary number of spacetime dimensions. Prom the result of that work, 
along with the previously obtained ODEs for tensor and vector perturbations,^ 
we have the full set of equations for perturbations of higher-dimensional maximally 
symmetric black holes in the form of a single self-adjoint 2nd-order ODE — which 
we refer to as the master equation — for each tensorial type of perturbation. These 
master equations, of course, can be applied to the perturbation analysis of higher- 
dimensional Schwarzschild black holes. Let us mention, furthermore, that these 
master equations have been derived in the framework of the gauge-invariant formal- 
ism for perturbations, and for this reason, they do not involve the problem of the 
choice of gauge. 

The purpose of this paper is to show that higher-dimensional Schwarzschild black 
holes are stable with respect to linear perturbations. As mentioned above, we have 
already obtained the master equation for each tensorial type of perturbation in the 
Schwarzschild background. Our main focus in this paper is therefore to prove that 
there exist no unstable solutions to the master equation with physically acceptable 
initial data. We also show that there exist no well-behaved static perturbations 
that are regular everywhere outside the event horizon. This is consistent with the 
uniqueness^ of a certain class of black holes in an asymptotically flat background. 
Although our main interest is in the standard higher-dimensional Schwarzschild black 
holes, that is, asymptotically fiat,*) spherically symmetric, static vacuum black holes 
in higher dimensions, we also investigate the stability of other maximally symmetric 
black holes with a non-vanishing cosmological constant with respect to tensor- and 
vector-type perturbations. We also give a stability criterion for generalised static 
black holes considered by Gibbons and Hartnoll.^ 

In the next section, we first give a general explanation of our background metrics, 
our formalism for gravitational perturbations, and the basic strategy of our stability 
analysis. At this stage, our arguments are not limited to the case of Schwarzschild 
black holes, but are equally applicable to any case of maximally symmetric black 
holes with a non-vanishing cosmological constant. Then, we perform the stability 
analysis for tensor, vector and scalar perturbations separately. In § |3J we exam- 
ine static scalar-type perturbations in the Schwarzschild background. Section 0] is 
devoted to summary and discussion. 

§2. Stability analysis 

As our background, we consider the (n + 2)-dimensional metric of the form 

dv 

ds 2 = -f{r)dt 2 + -—+r 2 dol, (2-1) 



*' It has been pointed out, however, that the concept of asymptotic flatness itself involves 
some subtle problems in higher dimensions, especially in the case of an odd number spacetime 
dimensions.^ In this paper, however, we do not discuss this issue and do not distinguish the cases 
of odd and even numbers of spacetime dimensions. 



4 



where 



/(r) = K-^L-\r\ (2-2) 



with constant parameters M, A and K = 0, ±1, and where cicr^ = ^{z)ijdz l dz^ 
denotes the metric of the n-dimensional maximally symmetric space K, n with sec- 
tional curvature K = 0, ±1.*) The constant parameter A is related to the (n + 2)- 
dimensional cosmological constant A as A = 2A/n(n + 1), and the parameter M 
defines the black hole mass. For example, in the case K = 1 and A = 0, this mass is 
given by 

nMA n 
8irG n+ 2 ' 

where A n = 2^ n+1 ^ 2 /r[{n + f )/2] is the area of a unit n-sphere, and G n +2 denotes 
the n + 2-dimensional Newton constant. We assume that M > to avoid the 
appearance of a naked singularity in our background. We are interested in the static 
region — called the Schwarzschild wedge — in which f(r) > 0, so that t is a time- like 
coordinate. Thus, when A > 0, the sectional curvature K must be positive, while 
when A < 0, K can be or ±1. The higher-dimensional Schwarzschild black hole 
corresponds to the case K = 1, M > and A = 0. 

To study perturbations in this background, it is convenient to introduce har- 
monic tensors in the n-dimensional symmetric space (/C n ,d<r^) and expand the per- 
turbations in terms of them, as we see below. The type of the harmonic tensor 
defines the type of the perturbation. Combining the expansion coefficients, one 
can construct the gauge- invariant variables in the two-dimensional orbit spacetime 
spanned by the coordinates y a = (t,r). The linearised Einstein equations then re- 
duce to a set of equations for the gauge-invariant quantities. It is found from the 
Einstein equations that there is a master scalar variable & for each tensorial type of 
perturbation. Consequently, the linearised Einstein equations reduce to a set of three 
master equations for the master variable. Thorough studies of the gauge-invariant 
formalism for perturbations and the derivation of the master equations are given in 
Refe.ES andHTJ). 

Because our background metric is static, with a Killing parameter t in the 
Schwarzschild wedge, the master equations can be cast into the form 

d 2 ( d 2 \ 

w*=[M- v r' (24) 

where r* = J drf~ l (r) and V is a smooth function of r*. Our task is then to show 
that the differential operator 

A =-£i +v < 2 ' 5 ' 

*' Accordingly, we also refer to the metric 12-11 1 as a maximally symmetric black hole One 
can also construct static, vacuum black hole metrics by simply replacing = ^ i jdz l dz : ' in 12-11 
by an arbitrary metric of an n-dimensional Einstein manifolo™*^**^* whose Ricci curvature is 
Rij — Kin — l)7ij. We briefly comment on such a generalised black hole case in §2] 
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is a positive self-adjoint operator in L 2 (r*, dr*), the Hilbert space of square integrable 
functions of r*, so that Eq. ()2-4j) possesses no negative mode solutions that are 
normalisable. In other words, the positivity of the self-adjoint operator ensures that, 
given well-behaved initial data — for instance, smooth data of compact support in 
r* — remains bounded for all time. A rigorous proof of the stability of the system 
described by Eq. (|2-4j) has been given, in the language of spectral theory, by WaldP^ 

It is known that, in the four-dimensional case, V — called the Regge- Wheeler 
potential for axial perturbations and the Zerilli potential for polar perturbations — 
is positive definite. Hence, the positivity of the self-adjoint operator immediately 
follows. However, in systems of greater than four dimensions, V is, as we will see, 
in general not positive definite. Thus it is not clear whether A can be extended to a 
positive self-adjoint operator. However, because V is bounded from below in most 
cases, in the domain Co°(r*) — the set of smooth functions with compact support in 
r* — A is a symmetric, semi-bounded operator and hence can be extended to a semi- 
bounded self-adjoint operator in such a way that the lower bound of the spectrum 
of the extension is the same as that of A. In fact, a natural extension, called the 
Friedrichs extension,^ is such a self-adjoint operator. Furthermore, when A > 0, 
there is a unique self-adjoint extension because in this case, the Schwarzschild wedge 
is globally hyperbolic, and the range of r* is complete, i.e., (—00,00). Therefore 
the time evolution of ^ — which is defined by a self-adjoint extension of A — must 
be uniquely determined by its initial data, by virtue of this global hyperbolicity. 
(If there were different self-adjoint extensions, then each of them would describe a 
different time evolution of & from the same initial data.) Therefore, if the symmetric 
operator A with domain Co°(r*) is shown to be positive definite, then the self-adjoint 
extension must also be positive definite. We show that this is indeed the case. 

When A < 0, the situation is quite different, because the range of r* is not com- 
plete, being (— 00, 0), and the Schwarzschild wedge is no longer globally hyperbolic. 
Nevertheless, there is a prescription's for defining the dynamics of in terms of 
a self-adjoint extension of A. We choose the self-adjoint extension that is obtained 
by imposing the Dirichlet boundary condition <P = at r* = as a natural require- 
ment for perturbations that are well-behaved (in r*). The other possible choices of 
the boundary conditions at r* = and the corresponding self-adjoint extensions in 
anti-de Sitter spacetime will be discussed elsewhere.^ 

In the subsequent sections, we examine tensor, vector and scalar perturbations 
separately, briefly reviewing the metric perturbations and the master equation. 

2.1. Tensor perturbations 

We begin by considering tensor perturbations, which are given by 



where is a function of y a = (t,r), and the harmonic tensors Tjj are defined as 
solutions to the eigenvalue problem on the n-sphere, i.e., 



Sgab = 0, Sg ai 



0, 5g l3 = 2r 2 H T T ij , 



(2-6) 



(A n + k&Tij = 



(2-7) 
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with 

Ti = 0, DjVi = . (2-8) 

Here, Dj is the covariant derivative with respect to the metric 7^ and A n = rf^DiDj. 
The eigenvalues are given by k T = l(l + n — l) — 2 for K = 1 and form a continuous 
set of non-negative values for K < 0. Note that in the expansion ()2-6|) . we have 
omitted the indices k labelling the harmonics and the summation over them. Note 
also that, in the four-dimensional spacetime (n = 2), there are no corresponding 
modes, as there are no harmonic tensors of this type on a 2-sphere. 

It can immediately be seen that the expansion coefficient Ht is itself gauge in- 
dependent, and it can be taken as the master variable for a tensor perturbation. The 
master equation follows from the vacuum Einstein equations [Eq. (5.5) in Ref. I26[) 
or, equivalently, Eq. (137) in Ref. l2*7|)] 

d ~t)^ =0, (2 ' 9) 

where <P = r n / 2 HT, and □ denotes the d'Alembertian operator in the two-dimensional 
orbit spacetime with the metric g a bdy a dy b = —f(r)dt 2 + / _1 (r)(ir 2 . 

The master equation (|2-9j) can be written in the form of Eq. (|2-4jl . with 



f 

V = V T = ^2 



(2-10) 



For the case of the standard Schwarzschild black hole (K = 1, A = 0), the asymptotic 
behaviour of Vp at large distances (r* — > 00) and near the horizon (r* — » —00) is 

C 

Vt ~ —j" for r* — > +00 , Vt ~ C_ exp (2Kr*) for r* — > —00 , (2-11) 

where k = (n — l)/2r^, with r^" 1 = 2M, and C± are positive constants. Typical 
profiles of the potential Vt are plotted in Figs. Q [21 and|31 We immediately see that 
the potential Vp is positive definite in the Schwarzschild wedge, irrespective of the 
signs of K and A [noting that, for K = 1, we have fcf. — (n — 2)K = (I — + n) > 0]. 
Thus, the operator A with the domain Cq ^*) is positive and symmetric (Hermitian) 
with respect to the inner product 

(#l,# 2 ) L2 = J du$\{r) <P 2 (r) . (2-12) 

Therefore, for the case r* G (—00, 00) (that is, the case in which A > 0), A is extended 
to a unique positive, self-adjoint operator in L 2 (r*). For the case A < 0, r* has the 
range (— 00, 0). For V = Vp, however, A is of the limit-point type*) at r* = 0, 

*' This type is characterised by the condition that for both ±i, one of the solutions to (— d 2 /dr^ + 
V)$ = ±i& is not square integrable at r„ =0. If A is of the limit-point type at both boundaries, 
then A — —d 2 /dr 2 + V is essentially self-adjoint. This can be seen by noting that Vt ~ n(n+2)/r 2 > 
3/(4rJ) near r, = oP 1 
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Fig. 2. The potential barriers Vt surrounding (n + 2)-dimensional Schwarzschild-de Sitter black 
holes for tensor perturbations of the mode I = 2. The horizontal axis is normalised with respect 
to the black hole horizon radius rn- 
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Fig. 3. The potential barriers Vt surrounding (n+2)-dimensional Schwarzschild-anti-de Sitter black 
holes for tensor perturbations of the mode I = 2. The horizontal axis is normalised with respect 
to the black hole horizon radius rn- 



and hence A is uniquely extended to a positive, self-adjoint operator. We therefore 
conclude that not only the standard Schwarzschild black hole (K = 1, A = 0) but 
all the possible maximally symmetric black holes in higher dimensions are stable 
with respect to tensor perturbations. This result agrees with that of the previous 
analysis P) 

2.2. Vector perturbations 

Next, consider vector perturbations, which are given by 

Sg ab = 0, Sg ai = rf a Yi, 6 9ij = -^-^H T D^ S) , (2-13) 

where the vector harmonics are introduced as the solutions of 

(A n + 4)Vi = 0, DiV = 0, (2-14) 

with eigenvalues k v that are given by k v = 1(1 + n — 1) — 1 with / > 1 for K = 1 
and form a continuous set of non-negative values for K < 0. 

A gauge-invariant variable F a can be constructed by combining the expansion 
coefficients f a and Ht- Substituting this into the Einstein equations, we can obtain 
the master scalar variable <P, in terms of which the gauge-invariant variable of vector 
perturbations for k v > (n — 1)K is given by 

F a = fa + ^~D a H T = r-^e ah D\r n l 2 ^) , (2-15) 
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where e a b denotes the anti-symmetric tensor on the two-dimensional orbit spacetime. 
Note that for K = 1 and I = 1, Vj becomes the Killing vector on the re-sphere. In 
this case, only perturbations corresponding to rotations of the background geometry 
satisfy the Einstein equations, as explained in Ref,l27|). Therefore we do not consider 
the I = 1 mode for K = 1 in what follows. 

The evolution part of the Einstein equations can then be reduced to [Eq. (5.14) 
in Ref. ESJ or, equivalently, Eq. (141) in Ref. EZJ] 

□ -^0 = 0, (2-16) 



which is identical to Eq. (|2-4|) with 



V = y v = ^ 



,2 / ,w . n(n + 2) n df 
k v -(n-l)K+^—f--r- 



(2-17) 



It can be shown that for the four-dimensional Schwarzschild case (n = 2, K = 1 
and A = 0), Vy coincides with the Regge- Wheeler potential. Also, it should be 
noted that the potential Vy falls off at large distances and near the horizon in the 
same manners as (j2-ll|) for K = 1 and A = 0. However, Vy is not positive definite 
for higher-dimensional Schwarzschild black holes with n > 2, K = 1 and A = 0. In 
particular, the constant C_ becomes negative when 2n > 21 — 1 + ^121(1 — 1) + 1 > 4 
(for / > 2), as seen in Fig.0J We can also see that the potential is not positive definite 
for the case of a non- vanishing cosmological constant in Figs. El and 

Now, to show the positivity of A, let us introduce the differential operator 

D = ^+S, (2-18) 

with S being some function of r. Then, we have 



where 









+ / 




boundary J 



(#, A$) L 2 = - $*D<P + dr*l\D$\ 2 + V\$r) , (2-19) 



V = V + - S 2 . (2-20) 
dr 

For <P £ C , ^°(r !t ), the boundary term vanishes. Inspecting the form of (|2-17|) . we find 



that the choice 



S=^ (2-21) 
2r v 1 



yields an effective potential that is manifestly positive: 

V = ^[k v -{n-l)K] >0. (2-22) 

Here, we note that for K = 1, k v — (n — l)K = (I — l)(Z + re) > 0. Hence, a symmetric 
operator A with domain C Y ^°(r^) is positive definite for any K and A (as long as the 
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Fig. 4. The potential barriers Vv surrounding (n + 2)-dimensional Schwarzschild black holes for 
vector perturbations of the mode I = 2. The horizontal axis is normalised with respect to the 
black hole horizon radius th- 

Schwarzschild wedge exits), so is the self-adjoint extension. For the case A < 0, the 
range of r* is incomplete. However, for the higher dimensional case (n > 3), V = Vy 
is shown to be in the limit-point case at r* = 0, and A is uniquely extended to a 
positive, self-adjoint operator. For the case of four spacetime dimensions (n = 2), 
imposing the Dirichlet boundary condition at r* = 0, as stated above, we can obtain 
a positive, self-adjoint extension of A. [Note that under the Dirichlet boundary 
condition, the boundary term in Eq. ()2-19|) does not contribute to the spectrum.] 
We have therefore confirmed that all the possible maximally symmetric black holes 
are stable with respect to vector perturbations. 

2.3. Scalar perturbations 

Scalar perturbations are given by 

figab = fab®, Sg ai = rf a Si, 5gij = 2r 2 (H L -yij§ + H T §ij) , (2-23) 

where the scalar harmonics §, the associated scalar harmonic vector §j, and the 
tensor S$j are defined by 

(A n + fc|)S = 0, Si = -— AS, S y = T 5 -A4- s + -7yS J (2-24) 

K§ Kg Tl 

with the eigenvalues kg given by kg = l(l+n — l) for K = 1 and forming a continuous 
set of non-negative values for K < 0. By definition, Sjj is a traceless tensor. Using 
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A,=0.05 
K=1, l=2 




Fig. 5. The potential barriers Vy surrounding (n + 2)-dimensional Schwarzschild-de Sitter black 
holes for vector perturbations of the mode / = 2. The horizontal axis is normalised with respect 
to the black hole horizon radius rn- 



the gauge-dependent vector defined on the orbit spacetime by 

r ( r 
ks V k S 

we can obtain the gauge-invariant variables for the scalar perturbation, 



(2-25) 



F = H L + -H T + -(D a r)X a , F ab = f ab + D a X b + D b X a . (2-26) 
n r 

Hereafter, we restrict our attention to modes for which / > 2, because the / = 1 modes 
do not correspond to any physical degrees of freedom, as discussed in Ref. I26|). 

We find, after a long calculation, that with the following choice of the gauge- 
invariant master variable ^, the linearised Einstein equations can be reduced to the 
form of Eq. (E3J) [or Eq. (3.11) in Ref.lSfi])]: 



nZ - r(X + Y) 



„n/2-l { m + n (n + l)x/2} ' 



where here and hereafter 



2M 



j , m = k s — nK , 



(2-27) 
(2-28) 

(2-29) 
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Fig. 6. The potential barriers Vv surrounding (n+2)-dimensional Schwarzschild-anti-de Sitter black 
holes for vector perturbations of the mode I = 2. The horizontal axis is normalised with respect 
to the black hole horizon radius rn- 

and X, Y and Z are denned by 

X = r n ~ 2 {Fl -2F) , Y = r n ~ 2 (F r r - 2F) , 8 t Z = -r n ~ 2 F t r . (2-30) 

The variables X, Y and Z can be expressed in terms of <P and its derivatives with 
respect to D a [see Eq. (3.12) and Eq. (3.10) in Ref.EEH]- 

The potential term in the master equation (|2-27|) is given by 



f(r)Q(r) 
16r 2 {m + n{n + l)x/2} z 



where 



Q(r) = - [n 3 (n + 2)(n + l) 2 x 2 - 12n 2 (n + l)(n - 2)mx + 4(n - 2)(n - 4)m 2 ] Ar 2 

+n 4 (n + l) 2 x 3 + n(n + 1) [4(2n 2 - 3n + 4)m + n(n - 2)(n - 4)(n + l)if] x 2 
-12n [(n - 4)m + n(n + l)(n - 2)] mx + 16m 3 + 4n(n + 2)m 2 . (2-32) 

The manners in which V$ falls off at large distances and near the horizon are the 
same as in the vector perturbation case. It can be shown that the potential Vs 
coincides with the Zerilli potential and is positive definite in the four-dimensional 
Schwarzschild case (n = 2, K = 1 and A = 0). However, for higher dimensions, 
the potential Vs is not positive definite in general, as shown in Fig. [7| (For further 
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Fig. 7. The potential barriers Vs surrounding (n + 2)-dimensional Schwarzschild black holes for 
scalar perturbations of the mode I = 2. The horizontal axis is normalised with respect to the 
black hole horizon radius r^. In 8-dimensional spacetimes, Vs is not positive definite. 



reference, we also plot the potential Vs for the Schwarzschild-de Sitter case in Fig. |H] 
and the Schwarzschild- anti-de Sitter case in Fig. 

From this point, we focus on the standard higher-dimensional Schwarzschild 
black hole case, i.e., that with A = and K = 1. Taking S to be 

S = £^, h = r n / 2+w {(1-1) (I + n) + n(n + 1W2} , (2-33) 
h dr 

we obtain as the effective potential defined in Eq. (|2-20|) the form 

V = f^)Q(r) (2 . 34) 

2r 2 {(/ - 1)(Z + n) + n(n + l)x/2} ' 1 ' 

where 

Q{r) = lx[ln{n + l)x + 2(1 - l){n 2 + n(3Z - 2) + (Z - l) 2 }] . (2-35) 

Clearly V > 0, and thus A with domain CQ°(r*) is strictly positive, so is the self- 
adjoint extension. The higher-dimensional Schwarzschild black holes are therefore 
stable with respect to scalar perturbations. This completes the stability analysis for 
higher-dimensional Schwarzschild black holes. 
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§3. Static perturbations 

In this section, we treat "static" scalar perturbations in the background of the 
Schwarzschild metric, with K = 1, A = and n > 2. Because scalar perturbations 
represent deviations from a spherical form, if there exist any regular static pertur- 
bations that fall off at large distances in such a manner that asymptotic flatness is 
maintained, it would be implied that there exist asymptotically flat, static vacuum 
black holes whose event horizons are topologically S n but are not spherical. We show 
that this is not the case. 

Regarding Y in Eq. (|2-3U|) as the master variable for a static perturbation, 
we have the following master equation for Y [with the replacement of r by x = 
2M/r n ~ l = (r h /r) n - 1 in Eq. (4.4) of Ref.HU)]: 

+ 2(n-2)x + 2 d + -(f-i)(Z + n) + ( w -2)x l y = Q 
dx 2 (n — l)x{x — 1) dx (n — 1) 2 j; 2 (1 — x) 



There are general solutions to Eq. (|.S-1|) that can be expressed in terms of hyper- 
geometric functions as 



15 




Fig. 9. The potential barriers Vs surrounding (n+2)-dimensional Schwarzschild-anti-de Sitter black 
holes for scalar perturbations of the mode I — 2. The horizontal axis is normalised with respect 
to the black hole horizon radius ru- When the number of spacetime dimensions is greater than 
4, Vs is not positive definite. 



+C2 py + v("- i+ ',- 2 <"- i ) + ',- 2 <" + '- i ) ;i y (3 .2) 

v r / \ n — 1 n — \ n — 1 / 

From this expression, we find that the original gauge-invariant variables F and F£ 
behave at the spatial infinity as pC\r l + qC2/r n+l ~ 1 , where p and q are non- vanishing 
constants that depend on the variables. Hence, in studying the question of whether 
regular, well-behaved static perturbations exist or not, we can choose the boundary 
condition C\ = 0, so that the gauge-invariant variable Y for I > 1 vanishes as r* — ► oo 
in the form 

Y - ° (^) • ( 3 - 3 ) 

Note that the I = modes preserve the spherical symmetry of the black hole, and 
hence need not be considered here. 

Now, we redefine the master variable as 



for which the asymptotic condition is written 

' 1 

r n/2+2-l 
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Then, we can show that Y obeys [see Eq. (4.8) in Eef. I2fi|) ] the equation 



■2 



iv- - (~+r ! v -i). ch„ 



with the potential [see Eq. (4.11) in Ref. 26)]: 



~ _ n(n - 2) + 41(1 + n - 1)(1 - x) + 2nx(l - x) + x 2 

Because this potential is positive definite, we can show that there exists no solution 
that satisfies the asymptotic condition (|3-6|) and is bounded at the horizon. First, 
from Eq. (|.3-6|) . there is a point r* = a at which Y is positive and dY/dr* = C < 0, 
by changing the sign of Y if necessary. Suppose that Y is not positive definite in 
r < a. Then, there must be a point r* = b < a such that Y is positive in the 
interval b < < a, and vanishes at r* = b. However, this cannot happen because 
d 2 Y/dr 2 > in b < r* < a and dY/dr* < at r* = a. Hence, Y" is always positive 
in r* < a. From this and the positivity of V, it follows that dY/dr* < C in r* < a. 
Hence, Y increases without bound as r* — > — oo. We therefore conclude that there 
exists no static perturbation that falls off at large distances and is regular everywhere 
outside the event horizon. 



§4. Summary and discussion 



We have shown that higher-dimensional Schwarzschild black holes are stable with 
respect to all tensorial types of linear perturbations, examining the master equations 
recently derived in the framework of the gauge-invariant formalism for black hole per- 
turbationsP^'^ We have confirmed that the 2nd-order elliptic differential operator 
A appearing in the master equation is a positive, self-adjoint operator in the Hilbert 
space L 2 , and we have thereby shown that each type of master equation possesses no 
negative-mode solutions that are normalisable. We have also proved that there exists 
no regular, well-behaved static scalar perturbation of Schwarzschild black holes. This 
result is consistent with the uniqueness of higher-dimensional spherically symmetric, 
static vacuum black holes. 

Our strategy explained in § [21 is not limited to the standard Schwarzschild black 
hole case. In fact, with the same method, we were also able to show the stability of 
maximally symmetric black holes for any K and A (and thus, including the higher- 
dimensional Schwarzschild-de Sitter and Schwarzschild-anti-de Sitter black holes) 
with respect to tensor and vector perturbations. Thus, to complete the stability 
analysis of maximally symmetric black holes, it is now necessary only to consider 
the case of scalar perturbations of non-asymptotically flat black holes. (Note that in 
the four-dimensional case, the stability of Schwarzschild-(anti-)de Sitter black holes 
has already been proved.^) 

The background metric (|2Tj) describes a generalised black hole when da^ is the 
metric of a generic Einstein manifold with Ricci curvature Rij = K(n — 1)7«. The 
stability of such a "black-Einstein" has recently been studied by Gibbons et a/P^^^ 
and HartnollP^' As a reason to concentrate on the analysis of tensor perturbations, 
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they argued that a generalised black hole should be stable with respect to vector and 
scalar perturbations because the standard higher-dimensional Schwarzschild black 
hole is believed to be stable with respect to these perturbations, at least in the case 
in which the event horizon takes the form of a compact Einstein manifold of positive 
curvature. Their argument is based on the following two points. First, because the 
Weyl tensor does not appear in the perturbation equations for the scalar and vector 
perturbations, the stability with respect to these modes is determined only by the 
spectrum of the Laplacian on the Einstein manifold. Second, the smallest eigenvalues 
of the Laplacian for these modes are always greater than or equal to those for the 
sphere, and hence the standard Schwarzschild black hole is the most unstable. Here, 
we regard the second smallest eigenvalue as the smallest one for the scalar mode, 
because the k = mode does not affect the stability analysis. Although they proved 
the second claim only for a vector mode, we can also prove it for a scalar mode with 
the help of the identity 

J dtTr&'Sij = ^~(k 2 -n) J da n SiS< > 0. (4-1) 

Furthermore, the first point implies that even for the generalised black hole case, the 
master equations we have examined here are still valid, if we replace the eigenvalues 
ky and kg of the Laplacian on K7 1 by the corresponding ones on the Einstein mani- 
fold.*) Our results in this paper thus place their analysis on a firmer basis, giving a 
proof for their argument summarised above for the case K = 1. 

We now discuss a few points regarding the stability of a generalised black hole 
with respect to tensor perturbations. In this case, the perturbation couples to the 
Weyl curvature. However, the master equation is still valid if we replace k 2 -, by 
\l — 2nK, with Xl being the eigenvalue of the Lichnerowicz operator on the Einstein 
manifold, as shown by Gibbons and Hartnoll. 11 ' Hence, from the positivity of Vr 
expressed in the form (|2-1U|) . we obtain a sufficient condition for the stability, 

Al > (3n - 2)K. (4-2) 



By rewriting Vr as 
V T = 

we obtain another sufficient condition 



n 2 - 10 + 8 n(n + 2) o n 2 M 



(4-3) 



. n 2 -10n + 8^ ,. 4 . 

X L > ~ A K, (4-4) 

for A < 0. This condition coincides with that derived by Gibbons and Hartnoll. 
Furthermore, if we apply the same trick used for vector perturbations to tensor 
perturbation by setting 

(4.5) 



The authors thank Sean Hartnoll for discussion of this pointP^ 
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(4-6) 



From this expression, we obtain the third sufficient condition, 



Al > 2{n-\)K. 



(4-7) 



For K = 1 and A < 0, the second condition l|4-4|l . is the strongest, while for K = 1 
and A > 0, the third condition, (|4-7[) . is the strongest. For K = —1, the first 
condition, (|4-2j) . is the strongest. 

Finally, we note that in the stability analysis carried out in the present paper, the 
supports of the initial data for the perturbations are, as one condition for being well- 
behaved, restricted to compact regions of the initial hypersurface inside the exterior 
Schwarzschild wedges: We have considered such initial data that vanish in an open 
neighbourhood of the event horizon, in particular, on the bifurcation n-sphere. In 
the four-dimensional case, this restriction on the initial data of the perturbations 
was removed by Kay and WaldP^ The key ingredients of Kay and Wald's theorem 
are the discrete isometry of extended Kruskal spacetime — which is implemented even 
in the higher-dimensional case (and also for cosmological horizons) — and the well- 
posed nature of the Cauchy problem for the wave equation — which also holds for 
the master equations. Furthermore, the forms in which the potentials V fall off near 
the horizon and in the asymptotic region are the same as in the four-dimensional 
case. Hence, there seems to be no dimensional obstruction in applying Kay and 
Wald's theorem to our present situation, which would imply the stronger claim that 
higher-dimensional Schwarzschild black holes are stable with respect to the linear 
perturbations whose initial data have compact support on a "Cauchy surface" in the 
extended Kruskal spacetime. 



The authors wish to thank Stefan Hollands and Bob Wald for many useful 
discussions and suggestions at the early stage of this work. We also wish to thank 
Gary Gibbons and Sean Hartnoll for discussions and comments, especially on the 
extension of the present analysis to the generalised black hole case. A.I. thanks the 
people of the Enrico Fermi Institute for their hospitality. This work was supported 
in part by the Japan Society for the Promotion of Science (AI), and by the JSPS 
grant no. 15540267 (HK). 



1) R. C. Myers and M. J. Perry, Annals of Phys. 172 (1986) 304 

2) E. Witten, Adv. Theor. Math. Phys. 2 (1998) 505 

3) S. Dimopoulos and G. Landsberg, Phys. Rev. Lett. 87 (2001) 161602 

4) S. B. Giddings and S. Thomas, Phys. Rev. D 65 (2002) 056010 

5) W. Israel, Phys. Rev. 164 (1967) 1776 

6) B. Carter, Phys. Rev. Lett. 26 (1971) 331 

7) D. C. Robinson, Phys. Rev. Lett. 34 (1975) 905 



Acknowledgements 



References 



19 



G. W. Gibbons, D. Ida and T. Shiromizu, Phys. Rev. Lett. 89 (2002) 041101 

G. W. Gibbons, D. Ida and T. Shiromizu, Phys. Rev. D 66 (2002) 044010 

R. Emparan and H. Reall, Phys. Rev. Lett. 88 (2002) 101101 

G. W. Gibbons and S. A. Hartnoll, Phys. Rev. D 66 (2002) 064024 

G.W. Gibbons, S.A. Hartnoll, and C.N. Pope, Phys. Rev. D 67 (2003) 084024 

G. Horowitz, and V. Hubeny, Phys. Rev. D62 (2000) 024027 

V. Cardoso, O. J. C. Dias and J. P. S. Lemos, Phys. Rev. D 67 (2003) 064026 
R. Konoplya, gr-qc/0303052 

C. Molina, gr-qc/0304053 

D. Ida, Y. Uch ida and Y. Mo risawa, Phys. Rev. D 67 (2003) 084019 
S. A. Hartnoll, hep-th/0305001 

F. R. Tangherlini, Nuovo Cimento 27 (1963) 636 

T. Regge and J. A. Wheeler, Phys. Rev. 108 (1957) 1063 

C. V. Vishveshwara, Phys. Rev. D 1 (1970) 2870 
R. Price, Phys. Rev. D5 (1972) 2419 

R. M. Wald, J. Math. Phys. 20 (1979) 1056 ;21 (1980) 218 

F. Zerilli, Phys. Rev. Lett. 24 (1970) 737 

F. Zerilli, Phys. Rev. D 2 (19 70) 2141 

H. Kodama and A. Ishibashi, hep-th/0305147 

H. Kodama, A. Ishibashi and O. Seto, Phys. Rev. D 62 (2000) 064022 
S. Hollands and A. Ishibashi, gr-qc/0304054 

D. Birmingham, Class. Quant. Grav. 16 (1999) 1197 

M. Reed and B. Simon, Fourier Analysis, Self-Adjointness, (Academic Press, New York, 
1975). 

A. Ishibashi and R .M. Wald, Class. Quantum Grav. 20 (2003) 3815 

A. Ishibashi and R. M. Wald, in preparation. 
S. A. Hartnoll, private communication. 

B. S. Kay and R. M. Wald, Class. Quantum Grav. 4 (1987) 893 



